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Solitons and fractional statistics 
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Solitons in the continuum limit of the Calogero model are derived and shown to correspond to one-particle 
excitations. The statistical mechanics of exclusion statistics particles is then formulated in terms of a priori 
probabilities and a path integral is thereoff constructed. 
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1. INTRODUCTION 

This talk consists, in fact, of two short stories, 
both connected to or motivated from the Calogero 
model. In the first, I derive an analytic expres- 
sion for a solitonic wave excitation in the contin- 
uum limit of the Calogero model, and show that 
it corresponds to one-particle excitations in the 
quantum description. Large-amplitude waves are 
also derived and correspond to a two-band quan- 
tum state. A conjecture for the chiral hamilto- 
nian of this problem is made. In the second, I 
formulate exclusion statistics in terms of micro- 
scopic occupation probabilities. This leads to a 
path-integral construction of the partition func- 
tion for such systems and can be used as a start- 
ing point for further generalizations. A bonus of 
this formalism is that the coefficients of the low- 
temperature expansion of the specific heat can be 
easily evaluated. Much of the material in this talk 
is taken from references . 

2. SOLITONS 

The Calogero-Sutherland-Moser model is ex- 
actly solvable in both the classical and the quan- 
tum regime This model is related to quan- 
tum spin chains with long range interactions be- 
tween the spins wave propagation in strat- 
ified fiuids 1^, random matrix theory and 
fractional statistics |q|. Remarkably, the quan- 



tum solution is much easier to interpret, exhibit- 
ing a straightforward analogy to the free fermion 
case. In a recent paper, Sutherland and Camp- 
bell examined the classical system in the thermo- 
dynamic limit and identified the excitations 
It was found that the classical system has soli- 
tons, corresponding to a single particle running 
through the rest of them, as well as small am- 
plitude waves (phonons), identified with holes. I 
shall derive here large amplitude wave and soli- 
ton solutions of the classical system in the con- 
tinuous limit, where the particles form a "fluid," 
and examine their correspondence to the quan- 
tum states. 

Consider a collection of particles of unit mass 
with the hamiltonian 



H 



N 



E 



(Xi XjY 



(1) 



We shall be interested in the limit N, L ^ oo 
with N/L fixed. In this limit, the system can be 
described in terms of a density field p{x) and a 
velocity field v{x). At equilibrium, the particles 
will form a regular lattice of spacing a and density 
Po = 1/a. The particle current is J = pv and by 
particle conservation 

p + dJ = p + d{pv) = (2) 

where d = d/dx. The kinetic energy of the sys- 
tem is 



'Part of this work was done while still at CERN 



K = J dx^pv^ 



(3) 
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We can formally solve eq. (2) for v to obtain 
V = —d^^p/p, and the expression for the kinetic 
energy becomes 

{d-^pf 



K ^ dx 



2p 



(4) 



This is exactly the kinetic term of the collective 
field hamiltonian description of a many-body sys- 
tem The potential energy can also be ex- 
pressed in terms of the density. The naive expres- 
sion, however is incorrect, the reason being that 
the interaction is singular at coincidence points, 
and thus a substantial part of the potential en- 
ergy comes from nearest neighbors. The correct 
expression requires a careful conversion of the dis- 
crete sum in terms of the continuous fields. Al- 
ternatively, we can simply take the classical limit 
{h —^ 0) of the quantum mechanical expression 
derived in the collective field formulation [TT| . 
The result is 



V 



dx 



^3 

6 ^ 



-pdp' 



9 {dpf 



2' » p 
where p stands for the Hilbert transform: 
1 



p= / dyP.P.- 



V 



p{y) 



(5) 



(6) 



The dynamics of the system can be found by 
varying the lagrangian L — K — V + fip with 
respect to p. The chemical potential p plays the 
role of a Lagrange multiplier ensuring that the 
total number of particles remains constant. The 
resulting equations of motion are 



+—P -9dp+-[-) 

as well as the continuity equation. By requiring 
that the static configuration u = 0, p = Po be 
a solution, we obtain the value of the chemical 
potential 

TT^ff 2 

^^-—Po 



(8) 



For a localized constant profile configuration, 
propagating at speed v, both p and v are func- 
tions of a; — vt only. From the continuity equation 
we have 



d{vp — yp) = and thus 



P- Po 



(9) 



In the above, the integration constant is fixed by 
the boundary condition that v at x ~^ ±oo, 
where p po- Similarly, the eq. of motion be- 
comes 



A solution of the above equation is 



Ap 



Pso\ = Po 1 + 



(11) 



{npox)^ + V? 

provided that v > v^. 

The above soliton carries particle number Q, 
momentum P and energy E, defined as the extra 
amount over the static solution po- We find 



Q 



dx (ps, 



Po) 



1 



P ^ dx Psol V = V 



(12) 



(13) 



E = Jdx [K{psoi) + y(psoi) - V{po)] = iv2(14) 

We observe that the net particle number carried 
by the soliton is 1, independently of its veloc- 
ity; its momentum and energy are also those of 
a free particle of unit mass moving at the soli- 
ton velocity v. Therefore, the soliton can be ex- 
actly identified with a particle excitation of the 
system. This is in agreement with exact results 
drawn from the quantum theory, where particle 
excitations always move faster than the sound ^ 
The form of the above soliton, however, is at odds 
with the results found in . We suspect that the 
source of the discrepancy is the truncation to a 
finite number of a;-derivatives of the form for the 
potential in |^ ; this turns the equation to a local 
one and gives the soliton an exponential decay, 
rather than the inverse-square decay of the non- 
local equation. 

A finite-amplitude periodic solution for the 
equations of motion is 



1 



1 



Pwavc Po+ , I , „ „ o 

A \^ VAM^ + l - XA COS ^ 

where 



-1 (15) 
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A is the amplitude of the wave and / its wave- 
length. The above equation is, therefore, the 
amplitude-dependent dispersion relation in this 
nonlinear system. Note that in the limit A oo 
the above equations reduce to the single soliton 
solution. 

In summary, we have found exact soliton and 
wave solutions for the CS system in the contin- 
uum limit. Certainly the above do not exhaust 
the list of solutions; the general motion of the sys- 
tem will be a nonlinear superposition of waves (or 
solitons). It is instructive to put the above solu- 
tions into correspondence with the quantum me- 
chanical states. Consider N particles in a space 
of length L. The ground state of the system 
consists of a "Luttinger sea" in the pseudomo- 
mentum, with spacing between adjacent particles 
equal to 2w£/L and "Fermi level" iriN/L, where 
g = h). At the limit H ^ 0, N,L ^ oo, 

N/L po, the ground state becomes a continu- 
ous filled band with Fermi level Pp = iTy^po. A 
small amplitude wave, corresponding to a hole, is 
a very small gap in the band. A soliton, cor- 
responding to a particle excitation, is a single 
particle peeled from the Fermi level and placed 
some distance above. The generic finite ampli- 
tude wave corresponds to a state with two con- 
tinuous filled bands, of widths Pi and P2 (with 
Pi+ P2 = ^Pp) and with a gap G between them. 
These are related to the wave parameters as 

A = (17) 

Such a state can be visualized as arising either by 
successively exciting single particles by the same 
constant momentum, until they form a continu- 
ous band, or by gradually augmenting the gap of 
a hole, until it becomes finite. This state can thus 
be thought of as either a coherent state of solitons, 
or as a coherent state of phonons, their nonlin- 
ear nature accounting for the change in profile as 
they accumulate. Indeed, the soliton itself can be 
thought of as a superposition of many phonons 
with very large wavenumber, and the phonon as 
a superposition of many solitons just above the 



Fermi level. For the finite N (finite L) system 
the distinction between the two is fuzzy and in 
principle only one kind of excitations need be 
considered as fundamental. Note, further, that 
quantum mechanically the holes behave as parti- 
cles with fractional statistics of order Ti/i (mean- 
ing that i/h of them put together would form a 
fermion). At the classical limit h ^ 0, thus, they 
become bosons, as they should be since phonons 
obey no exclusion principle. Particles, on the 
other hand, carry statistics of order i/h. Thus in 
the classical limit they become "superfermions" 
meaning that no two of them can occupy rel- 
atively nearby quantum states. This is consis- 
tent with the inverse square repulsion between the 
classical particles. It should be noted that a sim- 
ilar soliton solution exists in the fluid description 
of this problem keeping the h corrections [Q , and 
it would be intersting to flnd its interpretation. 

We conclude by noting that the quantum me- 
chanical problem separates into two noninteract- 
ing chiral sectors, having to do with excitations 
near either end of the Luttinger sea. (The two 
sectors mix nonperturbatively when a number of 
particles of order N is excited, depleting the sea.) 
Therefore, the equation governing the continuum 
system should also decompose into two nonmix- 
ing, first-order in time equations, one for each 
sector. For the corresponding equation for free 
fermions this is indeed the case jlj]. In fact, from 
the collective field description of the system when 
only one chiral sector is present, we deduce that 
the chiral equations are exactly of the Benjamin- 
Ono type |l^Jl^. The exact field combinations 
in terms of which this decomposition would be 
achieved, however, are not known and constitute 
an open problem. 

3. EXCLUSION STATISTICS 

Statistics is an inherently quantum mechani- 
cal property of identical particles which, as the 
name suggests, modifies the statistical mechani- 
cal properties of large collections of such particles. 
Motivated by properties of the Haldane-Shastry 
model (a lattice version of the Calogero model), 
Haldane defined a generalized exclusion statistics 
through the reduction of the Hilbert space of ad- 
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ditional particles in a system due to the ones al- 
ready present in the system |15|. He proposed 
then the definition 
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Ad 

"an 



(19) 



where N is the number of particles in the sys- 
tem, d is the dimensionality of the single-particle 
Hilbert space, obtained by holding the quantum 
numbers of — 1 particles fixed, and Ad and AA^ 
are their variation keeping the size and boundary 
conditions of the system fixed. 5 = corresponds 
to bosons (no exclusion) while g = 1 corresponds 
to fermions, excluding a single state for the re- 
maining particles, the one they occupy. 

On the basis of the above, Haldane proposed 
the combinatorial formula for the number of 
many-body states of N particles occupying a 
group of K states 



M = 



[K-{g-l){N-l)]l 
Nl[K - g{N - 1) - l]l 



(20) 



Starting from this, the thermodynamical proper- 
ties of particles with exclusion statistics can be 
derived and this system has received a 

lot of recent attention Jl9|"p^ . 

It is obvious from that exclusions statistics 
makes sense only in a statistical sence, since Ad 
and M can become fractional for AA^ = 1 or 
= 1. It is, nevertheless, useful to attempt a 
microscopic realization and interpretation of frac- 
tional exclusion statistics, and see what it implies 
for the one-particle states. Such a description has 
the obvious advantage of being generalizable to 
interacting particles, for which the notion of d 
becomes hard to define. 

The starting point will be the grand parti- 
tion function for exclusion-statistics particles in 
K states 



Z{K) = AI{K,N)a 



.N 



(21) 



N=0 



where we put x — exp(/i — e)/A;T with fi the chem- 
ical potential and assumed that all K states are 
at the same energy e. In the statistical limit of 
large K, Z{K) should be extensive. This intro- 
duces, then, the notion of a microscopic descrip- 
tion of the system in which the above Z is the 



K-ih power of a single-state partition function. 
Each single level can be occupied by any number 
of particles, but with an a priori probability P„ 
for each occupancy n independent of the temper- 
ature. We thus demand 

Z{K, x) = Pn{K)xA = M{K, N)x^{22) 

\ n J N=0 

for all X. The above probabilities must, in gen- 
eral, depend on K. If, however, Pn{K) assume 
some (finite) asymptotic values as K goes to in- 
finity (as they should for an extensive Z), then the 
above microscopic partition function becomes an 
accurate description in the statistical limit. This 
indeed happens, and we obtain 



(23) 



The most obvious feature of the above expres- 
sions is that, unless g = 0,1, they always become 
negative for some values of n. Therefore, their in- 
terpretation as probabilities is problematic. This 
is an inherent problem of fractional g-statistics 
which cannot be rectified by, e.g., truncating 
M {K, N) to zero for A^ > K/g. The description 
of the statistical system in terms of effective nega- 
tive microscopic probabilities is, nevertheless, ac- 
curate and useful. Note, also, that the above P„ 
never truncate to zero for n above some maximal 
value (unless g ^ I), unlike parafermions. 

From the above, the single-level partition func- 
tion Z(x) = Z can be shown to satisfy 



Z9 _ ^9-1 



(24) 



The average occupation number n is expressed as 
1 



n = — xdxZ — xdxW 
Z 



(25) 



where W = \nZ is the thermodynamic potential 
(over — kT). It can be shown that the above imply 



(1 - gny[l -{g- l)n]^-9 = fix^i 



(26) 



in accordance with the result of [[16[-|18[ . 

From the above we can easily calculate the low- 
temperature expansion coefficients of the energy 
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(or specific heat) in terms of tlie temperature. Ex- 
cept for factors depending on the dimensionahty 
of space, these are Q 

C„= / — (lna;)"[(-l)"n(^)+fi(x-i)--](27) 
Jo X g 

where 1/g is the saturation density for n at zero 
temperature and e < fj,. The trick is to use (25) 
and change variable of integration from x to W. 
Since by (24) x is an exphcit expression of W, we 
obtain the exphcit integral 

C„ = (-l)"y {[5W^ + ln(l-e-^)r 

° ~igWr}dW (28) 



The first few coefficients are 

Co^O, Ci-^, C2--C(3)(l + .g), etc. (29) 



We see that Co ~ for all g, as conjectured 
in This expresses the fact that the ground 

state of the many-body g-on system is nondegen- 
erate, which is expected to be a generic feature 
of particle systems. Ci is also independent of g. 
This result has first been conjectured by S. Isakov 
and then numerically verified by J. Myrheim. 
The first analytical derivation was provided by 
D. Arovas [Q. In fact, from (28) we see that 
Cn is a polynomial in g of degree n — 1, its fc-th 
coefficient being 



W^'=[ln(l-e-'^)]"-'=dW^(30) 



Through a change of variables we can show that 
the above coefficients obey 



C„ 



Cn 



(31) 



so only half of the coefficients need be calculated. 

The above relation points to a mapping be- 
tween the small-g and large-g regions. Indeed, 
this is related to a duality relation in these sys- 
tems. It is easy to see that Z satisfies: 



Z-\g,x-3) + Z-\g-\x) = l 



(32) 



From the above relation, the duality relation for 
the density is recovered flS] , ^ 

(33) 



gnig,x)+g-^n{g~\x-'/3-) = 1 



We regard the formula for Z (32) as more fun- 
damental since it seems to be more generic. For 
instance, parafermions of order p = l/g are de- 
fined such that at most p particles can be put per 
state with probabilities 1. Thus 



^par 1 3^ ■ 



1 - a;P+i 
l~x 



(34) 



from which we can write the generalized 
parafermionic partition function Zparig^x) by 
simply putting p = l/g above. It can be seen 
that Zpar also satisfies (32) but not (33). 

The free energy W can be expressed as a power 
series in x 



w^y^x- 

^ n 

n—l 



(35) 



in terms of the "connected" weights 

wi = PI, W2 - 2P2-P1', W3 = 3P3-3PiP2+Pf(36) 



etc. We find for w„: 



w„ 



n 



gn 
m 



(37) 



These are remarkably similar to P„ (except for 
the range of m) . Notice that the w„ are not prob- 
abilities, but rather cluster coefficients. In fact, 
Wn = 1 for bosons and w„ = (— )"~"'^ for fermions. 
Also, W2 = I - 2g 

From the above expressions for Wn we can find 
a path integral representation for the partition 
function of exclusion particles in an arbitrary ex- 
ternal potential. We start from the usual eu- 
clidean path integral with periodic time /3 for N 
particles with action the sum of N one-particle 
actions, and sum over all particle numbers N with 
appropriate chemical potential weights. Since the 
particles are identical, we must also sum over 
paths where particles have exchanged final po- 
sitions, with weights equal to the inverse symme- 
try factors of the permutation to avoid overcount- 
ing (compare with Feynman diagrams). Thus the 
path integral for each N decomposes into sectors 
labeled by the elements of the permutation group 
Perm(N). By the usual argument, the free en- 
ergy will be given by the sum of all connected 
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path integrals. It is obvious that these are the 
ones where the final positions of the particles are 
a cyclic permutation of the original ones (since 
these are the only elements of Perm{N) that can- 
not be written as a product of commuting ele- 
ments). These have a symmetry factor of 1/7V 
corresponding to cyclic relabelings of particle co- 
ordinates (compare with the factors of l/n in- 
cluded in (35)). They really correspond to one 
particle wrapping N times around euclidean time. 
Thus, if we weight these configurations with the 
extra factors wj^ , as we have the right to do since 
they belong to topologically distinct sectors, we 
will reproduce the free energy of a distribution of 
g-statistics particles on the energy levels of the 
one-body problem, that is 



is due to the above projection, which must be in- 
serted in the (unitary) many-flavor path integral. 

The possibility to define statistics through the 
choice of the coefficients w„ suggests other pos- 
sible generalizations. Perhaps the simplest one is 
to choose 



(40) 



that is, one factor of —a for each unavoidable 
particle crossing. This leads to the statistical dis- 
tribution for the average occupation number n 



1 



(41) 



oo . N 

N=l •' n=l 



'SE[Xr.it„)] 



where Se is the one-particle euclidean action and 
the paths obey the boundary conditions x„ (/3) = 
a::„+i(0), xnIP) = a;o(0). {x can be in arbi- 
trary dimensions.) The partition function will be 
the path integral over all (connected and discon- 
nected) paths, with appropriate symmetry factors 
and a factor of w„ for each connected rt-particle 
component. 

It is clear that the above path integral is not 
unitary, since the weights Wn are not phases, nor 
does it respect cluster decomposition, since the 
Wn do not provide true representations of the per- 
mutation group (unlike the g = 0,1 cases). This 
is again a manifestation of the non-microscopic 
nature of exclusion statistics. It does make sense, 
nevertheless, at the statistical limit. 

The above path-integral realization can be ex- 
tended to other statistics. E.g., for parafermions 
of order p (with p integer) the corresponding 
weigts Wn are 



which is the simplest imaginable generalization of 
the Fermi and Bose distribution, and was exten- 
sively analyzed in |24[ . The combinatorial for- 
(38) mula for putting N particles in K states for the 
above a-statistics is 

^ K{K - a){K -2a)---{K-{N- l)a) 

m ^ ' 

This can be thought as a different realization of 
the exclusion statistics idea: the first particle put 
in the system has K states to choose, the next has 
K — a due to the presence of the previous one an 
so on, and dividing by A^! avoids overcounting. 
Fcrmions and bosons correspond to a = 1 and 
a = — 1 respectively, while a — Q corresponds 
to Boltzmann statistics (as is also clear from the 
path integral, in which no configurations where 
particles have exchanged positions are allowed, 
but factors of 1/iV! are still included). The cor- 
responding single-level probabilities are 



Pn = 



n— 1 ^ 



-|- TO 



(43) 



Wn 



-p for n = 0mod(j3+l) , 1 otherwise. (39) 



This representation is more economical than the 
one calling for p distinct flavors of fermions 
and projecting into invariant states (singlets) un- 
der Perm(p). The origin of the apparent non- 
unitarity and breakdown of cluster decomposition 
in the above integral for parafermions is clear: it 



For a = 1/p with p integer (a fraction of a 
fermion), the above probabilities are all positive 
for n up to p and vanish beyond that. For a < 
all probabilities are positive and nonzero. Thus, 
the above system has a bosonic (a < 0) and a 
fermionic (a > 0) sector, with Boltzmann statis- 
tics as the separator. It is a plausible alterna- 
tive definition of exclusion statistics, due to (42), 
and has many appealing features, not shared by 
the standard (Haldane) exclusion statistics, such 
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as positive probabilities, a maximum single-level 
occupancy in accordance with the fraction of a 
fermion that a represents, and analytic expres- 
sions for all thermodynamic quantities. It would 
be interesting to find a physical system in which 
these statistics are realized. 

We conclude by pointing out that, once we have 
the path integral we can easily extend the notion 
of exclusion statistics to interacting particles: we 
simply replace the action ^E[xn] by the full 
interacting iV-particle action, thus circumventing 
all difficulties with combinatorial formulae. In 
the interacting case one has to work with the full 
partition function, rather than the free energy, 
since topologically disconnected diagrams are still 
dynamically connected through the interactions 
and do not factorize. Applications of the above 
on physical systems would be welcome. 
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